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ABSTRACT 

Some  problems  in  the  optimum  control  of  a  linear  dynamic 
system  are  investigated,  particularly  the  problem  of  deter- 
mining the  minimum  time  required  to  drive  a  linear,  constant- 
coefficient  dynamic  system  from  an  initial  state  to  a  specified 
terminal  state  with  a  limited  power  scurce.   The  important 
feature  of  the  paper  is  that  an  elementary  method  of  solution 
fcr  this  problem  is  given.   It  is  a  method  of  successive 
approximations  based  on  the  adjoint  system  of  differential 
equations  in  a  way  similar  to  that  which  bliss  used  in  cal- 
culating differentials  in  Ballistics.   A  program  is  given  for 
solving  the  minimum- time  problem  on  a  digital  computer;   an 
elementary  proof  is  given  that  if  the  routine  converges  then 
the  solution  thus  found  yields  the  desired  minimum  time. 

The  forcing  vector  is  assumed  to  have  bounded  components, 
representing  limited  driving  power  or  voltage.   This  forcing 
vector  is  to  be  selected  as  a  function  of  time  sc  as  to  drive 
the  dynamic  system  from  one  state  to  another  in  minimum  time. 
In  problems  where  certain  control  variables  are  bounded,  it 
is  well  known  that  the  control  function  assumes  at  all  times 
either  its  maximum  or  its  minimum  value,  changing  d is continuous- 
ly from  one  to  the  other.   This  is  generally  the  case  here. 
The  problems  v/ith  this  type  of  control  are  frequently  referred 
to  as  Bang -Bang  Control  problems  [l] ,[2] ,[3]2. 

1.  A  system  of  differential  equations  will  be  considered  as 
representing  a  dynamic  system. 

2.  References  to  the  bibliography  will  be  denoted  by  numbers 
in  brackets  [   ] . 
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The  solution  cf  the  minimum-time  problem  for  an  N-th 
order  linear,  constant-coefficient  system, in  which  the  matrix 
of  coefficients  has  real  and  distinct  eigenvalues,  is  given 
in  matrix  notation,  including  the  set  of  differential-correction 
equations  for  determining  the  adjoint  vectcr  and  minimum  time 
associated  with  the  desired  terminal  state.   The  solution  of 
the  adjoint  system  is  given  in  terms  of  exponential  functions 
of  the  eigenvalues  of  the  matrix  of  coefficients  end  the 
corresponding  eigenvectors  (see  Appendix  I),  the  components  of 
v/hich  are  in  the  same  ratio  as  the  cofactcrs  of  the  elements 
cf  a  column  cf  the  corresponding  characteristic  matrix  of  the 
original  system  (see  Appendix  II). 

The  general  case  of  the  second -order  minimum-time  problem 
was  programmed  and  run  on  the  CDC-1604  computer.  Because  the 
convergence  of  the  differential-correction  equations  may  be 
critical,  especially  for  stable  systems,  a  search  routine  was 
incorporated  into  the  program  to  search  the  input-parameter 
space  fcr  "useable"  starting  values,   A  second  supplementary 
program  provides  for  the  plotting  of  minimum-time  trajectories 
with  initial  state  X(0)  =  0  fcr  any  second-order  system. 

The  general  third -order  minimum-time  problem  was  also 
programmed  and  run  en  the  CDC-1604,  using  the  general  matrix 
techniques  described  for  the  N-th  order  system.   Extension  of 
the  program  tc  N-th  order  could  be  accomplished  directly. 

1.   A  minimum-time  trajectory  is  the  locus  of  points  in  the 
state  space  through  which  the  dynamic  system  passes  as  it  is 
driven  from  an  initial  state  to  a  terminal  state  in  mini  num 
time. 
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A  plot  cf  the  minimum-time  trajectories  with  initial 
state  X(0)  =  0  for  a  particular  stable  system  is  included 
in  illustration  20   This  illustration  also  contains  the 

switching  curves  associated  with  this  initial  state  and 

2 
the  region  of  static  stability,,   Switching  curves  do  not 

seem  to  have  a  very  general  significance,  i.e.,  they  are 

not  implicit  in  the  differential  equations  but  depend  on 

the  initial  conditions „   It  is  to  be  noted  however  that 

for  a  second-order  dynamic  system  and  those  minimum-time 

trajectories  with  two  sv/itching  times,  the  time  between  the 

switching  times  is  an  invariant  of  the  dynamic  system. 

The  method  of  differential  corrections  for  the  adjoint 
system  was  first  developed  by  Bliss  for  making  differential 
corrections  in  Ballistics  [8]0   The  set  of  differential- 
correction  equations  developed  in  this  paper  provide  for 
the  correcting  of  the  input  parameters  to  obtain  the  adjoint 
vector  and  minimum  time  associated  with  specified  initial 
and  terminal  states 0   None  of  the  papers  reviewed  on  the 
Bang-Bang  Control  Problem  provide  such  a  system, 

I  would  like  to  express  my  gratitude  for  the  encourage- 
ment, guidance,  and  inspiration  which  Professor  Faulkner  has 
provided,  and  without  which  this  paper  would  not  have  been 
pcssible8 

1.   A  switching  curve  (surface)  is  the  locus  of  points  at 
which  one  of  the  components  of  the  optimum  forcing  function 
changes  sign0 

20  #A  point  of  static  stability  of  the  system  represented 
by  X  =  AX  +  F  is  a  point  such  that  an  admissible  force  can 
be  chosen  to  constrain  the  system  to  that  point.   That  is,  it 
is  a  point  where  \£ia±-xA    <  1,  (i=l,2,. „ . ,n) ,  and  x±   =  0, 
(i=l,2, o . . ,n)  for  some  admissible  force  F   (see  section  4) 0 
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1.  Introduction. 

Let  us  consider  a  linear     tic  system  represented  by 
the  system  of  differential  equations  in  matrix  form 

(1)  X(t)   =  A  X(t)   +  F(t) 
or  componentwise  as 

(2)  x.(t)  =  £  a^x^t)  +  f.(t)      (i  =  l,2,...,n). 

k 

The  matrix  A  is  an  n  x  n  matrix  with  constant  elements  a., 
and  with  real,  distinct  eigenvalues e 

X(t)  is  an  n  x  1  matrix  representing  a  point  or  a  state 
of  the  dynamic  system.   Kc  distinction  will  be  made  generally 
between   an  n  x  1  (column)  matrix  and  a  vector.   If  a  quantity 
is  to  be  considered  as  a  vector,  as  in  a  scalar  product,  a 
bar  (— )  will  be  placed  over  it. 

The  dynamic  system  is  controlled  by  a  limited  power  source, 
i.e.,  the  components  f=(t)  of  the  vector  forcing  (control) 
function  F(t)  are  restricted  by  the  bounds 

(3)  |f::(t)|   <  L±  (i  =  l,2,...,n) 
where  the  L.  are  preassigned  constants.   The  L.  may  all  be 
assumed,  to  be  1  (one)  without  loss  of  generality,  and  this 
v/ill  be  the  case  henceforth . 

The  principal  problem  under  investigation  is  the  deter- 
aticn  of  the  minimum  time  required  to  drive  the  dynamic 
system  from  an  initial  state  X(0)  to  a  specified  terminal 
state  Xf,  and  the  determination  of  the  control  force  F(t) 
which  achieves  this  minimum  time  and  the  corresponding  path 
or  trajectory,, 


The  solution  is  effected  by  using  the  adjoint  system 
of  differential  equations.   Certain  properties  of  the  adjoint 
system  make  this  method  well  suited  for  the  solution  of  this 
problem. 

The  adjoint  system  of  differential  equations  is  a 
homogeneous  system;  the  system  adjoint  to  (1)  and  (2)  is 
formally  defined  by  Bliss  [9]  as 

(4)   U(t)   =   -  AT  U(t) 
where  the  matrix  A  need  not  have  constant  elements.   The 
superscript  T  denotes  the  transpose  of  a  matrix. 

Written  componentwise,  the  system  adjoint  to  (1)  and  (2) 
is 

(b)   u. (t)   =  -  Z   a,  ,u,  (t)  (i  =  l,2,..,,n). 

i  k  Kl  k 

Whenever  the  elements  of  A  are  constant  the  homogeneous 

adjoint  system  can  be  solved  directly  in  terms  of  exponential 

T 
functions  of  the  eigenvalues  of  the  matrix  of  coefficients, -A  , 

and  the  corresponding  eigenvectors. 

The  fundamental  adjoint  formula  [9]  of  Bliss^ 
t2 

allows  us,  by  proper  selection  of  the  adjoint  solutions 
(choosing  those  solutions  of  the  adjoint  system  which  meet 
specified  end  conditions),  to  express  the  terminal  values 
of  the  solutions  of  the  original  set  of  inhomcgeneous  equations 
in  terms  of  integrals  of  the  adjoint  solutions  and  the  forcing 
function  of  the  dynamic  system. 

1.  See  Appendix  I. 

2.  See  section  2  for  derivation  and  proof, 
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(6)    [  I   u.x.  ]t*    =  f  2  (  I  u.f.)  dt 


The  order  cf  development  in  this  paper  will  be  as 
follows: 

First,  the  adjoint  system  cf  differential  equations  is 
defined,  and  the  relation  is  established  between  solutions 
cf  the  original  system  of  equations  and  its  adjoint  through 
the  fundamental  adjoint  formula,  a  form  cf  Green's  Formula[ll]. 

Second,  the  minimum-time  problem  is  solved  by  applying 
the  forcing  vector  such  that  its  projection  on  a  prcperly- 
.cen  adjoint  vector  is  at  all  times  a  maximum.   This 
perty  is  also  known  as  the  Maximum  Principle  of  Pontryagin. 
The  solution  requires  the  solution  cf  a  two -point  boundary- 
value  problem.   The  method  employed  is  a  method  of  successive 
pprcximations,  using  the  differential  procedures  which  Bliss 
introduced  in  Ballistics,  based  on  the  adjoint  system  [8]0 
Further,  an  elementary  prccf  is  given  that  a  solution  so 
obtained  furnishes  the  minimum  time. 

Third,  an  example  is  given  of  the  solution  of  a  specific 
second-order  system  through  the  use  of  the  adjoint  system  cf 
differential  equations. 

Fourth,  the  solution  of  the  general  N-th  order  problem 
is  derived  in  matrix  notation  for  the  case  cf  a  linear,  constant- 
coefficient  system  of  differential  equations  where  the  matrix 
of  coefficients  has  real  and  distinct  eigenvalues.   This 
includes  the  differential-correction  equations  for  determining 
the  adjoint  vector  and  the  minimum  time  associated  with  a 
specified  terminal  state 0 


Fifth,  the  techniques  described  abcve  are  applied  in 
the  programming  en  a  digital  computer  of  the  solution  of  the 

i  rrum-time  problem  for  general  second-  and  third -order 
systems. 


The  Adjoint  System  cf  Differential  Eq       ,  Green's  Formul 
Let  is  consider  a  syste   cf  Line,  i   ifferential  equations 
j  tier:  ; o   mt  1"  i  .  foi 

(1)  X  =  A  X  +  F 
or  componentwise  ■ 

(2)  Xj  =  T.     -ik^,  +  f3  (i  =  l,2,...,n) 

K. 

ire  the  elements  a.,  of  mat  i;  A  are  constants  (or  continuous 

J_lv 

funetic  is  cf  t),  and  ..'here  f.  ore  functions  cf  t  [9]. 
Definition ; 

The  system  cf  differential  equations  adjoint  tc  (1)  and  (2) 
is  defined  [9]  as 

(3)  V      =      -  Ax  U 
cr  cc        i  se  as 

Tl    .  .  (7.   :  .ontgl  Adjoint  Formula)  [9]: 

If  f(2)   x.   =  I  a.kxk  +  f .  (i  =  l,2,...,n) 


a. 


>(4)   v.±     =  -  £  c,  ..;  (i  =  l,2,...,n) 

here  (4)  is  the  adjoint  of  (2), 

then  (5)  -£(£ux)  =   £V   . 

i  i 

t       t  held  for  every  solution  of  (2)  one  (4). 
Proof: 


C  i  .  -'  i  r  the  Lagrange  Identity  [H] 

(6)   d  /  v      N  v*  , 

dt  t  L  ^^M)  -  L     (u.x,  +  u.x, 

i  i    x  x    -1-  j- 


) 


Substituting  (2)  and  (4)  into  (6)  we  get 

ilk  i 

Expanding  (7)  we  get 

^  -nkTu.x.)  =  F I!  a.,u.x,  +  5  u.f.  -  T  T   a,  .u,.x.  , 

l        l  k  i        i  k 

and,  interchanging  the  meaning  cf  i  and  k  in  the  last  term, 

^t(  Z  u.x. )  =  £  i]a.,  u.x,  +  £  u.f .  -  T.  T  a.,  tux,  , 

d  t  v  i  i    .  r  i/,  i  k    T*  l  i    i  t  ik  i  k  ' 
i         i  k  l         k  l 

which  upon  interchanging  the  order  of  summation  in  the  last 

term  becomes 

(10)  4([ii.x.)   =  TT   a.,u.x,  +[u.f.  -  [[a.u.x,  . 
dt  V  i  l     *?*  t     lk  i  k   4-  i  i    V  ,   lk  l  k 

l         l  k  i         l  k 

Kence  we  get  the  final  result,  the  Fundamental  Adjoint  Formula 
cf  bliss  [9]: 

(5)  ;4(  Z   u.x.)  -  j>.f,    • 

o  t  7   i  i '     v  i  i 
l  l 

Alternatively,  the  Fundamental  Adjoint  Formula  can  be 

written  in  integral  form  as 

til)   [  £  u.x.  ]t*   =   f  '*    (  £  u.f.)  dt  . 
i  J  j_     i 

This  form  is  cften  referred  to  as  Green's  Formula  [11], 


3.  Outline  cf  Solution,  Preliminary  Formulas. 

Let  us  now  consider  the  system  cf  differential  equations 
discussed  in  section  1 

(1)   X(t)   =  A  X(t)   +  F(t)  . 
In  this  paper  the  dimension  of  F(t)  is  the  same  es   that  of 
X.(t),  though  this  is  irrelevant  to  the  method  cf  solution. 

Given  an  initial  state  X(0)  of  the  dynamic  system  and 
another  state  X~  (referred  to  as  the  terminal  state),  we 
desire  tc  determine  how  to  drive  the  system  from  X(0)  tc  X~ 
in  minimum  time  by  prcper  application  of  the  control  fcrce 
F(t). 

The  solution  cf  the  minimum-time  problem  consists 
basically  of  the  following  steps: 

a.  Equations  are  obtained  from  Green's  Formula   (2.11) 

tc  express  the  terminal  values  of  the  system  state  variables 
in  terms  of  the  initial  state  values,  the  time,  and  an 
integral  involving  the  fencing  function  F(t)  and  the  sclutions 
cf   the  adjoint  system „ 

b.  It  is  established  next  that  if  a  solution  is  obtained 
which  maximizes  a  certain  integral  and  satisfies  the  desired 
end    conditions  it  is  the  desired  solution. 

c.  Only  solutions  are  sought  which  maximize  an  integral 
of  this  form.   However,  there  are  seme  parameters  in  the 
integral,  and  these  must  be  determined  so  as  tc  yield  the 


1.   Equations  from  other  sections  will  be  denoted  by  decimal 
numbers;  the  number  preceding  the  decimal  point  indicates 
the  sect j on,  and  the  number  fc Hewing  indicates  the  equation 
number  in  that  section. 


desired  end  conditions.   The  problem  reduces  tc  the  solution 
cf  a  two-point  boundary- value  problem. 

d.   This  problem  is  solver;  by  a  method  of  successive 
approximations.   The.  corrections  tc  the  input  parameters 
are  obtained  by  a  method  based  on  Bliss'  method  of  computing 
differentials  in  Ballistics „ 

DeSoer  [2]  proves  the  existence  of  a  minimum-time  path 
for  stable  systems. 

The  solution  of  the  adjoint  system  obtained  in  the 
process  defines  the  functional  for  the  Maximum  Principle  of 
Pontryagin  [7].   The  integrand  cf  the  functional  is  the 
scalar  product  cf  the  forcing  vector  and  a  vectcr  defined 
by  the  solution  of  the  adjoint  system  of  differential  equations 

We  therefore  have  under  investigation  the.  dynamic  system 
governed  by  the  system  cf  differential  equations 

(1)  X  =  A  X  +  F  . 

The  system  adjoint  tc  (1)  can  be  written  directly  as 

(2)  U  =  -  A1  U0 

From  the  Fundamental  Adjoint  Formula  (Green's  Formula) 
we  get,  cheesing  tj  =  0,  t2  =  T, 

T     rT 

(3)  [  £  u.x.  ]Q  =  /    (  £  uifi)  dt  . 

1  o     1 

Let  us  now  examine  the  solutions  of  the  homogeneous 
adjoint  system.   From  Appendix  I  we  see  that  (for  n  distinct 
eigenvalues  -m.  of  matrix  -A  )  there  are  n  distinct  solutions 
L"J  cf  the  adjoint  system  cf  the  form 

(4)  U^  =  e'm^   Cj  (j  =  1,2,. ..,n) 
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where  CJ  is  the  eigenvector  [c,  .  cQ.  ....  c  .]x  as sedated 

with  eigenvalue  -m . .  U^  =  [uj  u£  .....  u^  ]  . 

Since  the  solutions  U^  are  independent  solutions  of 
the  adjoint  system,  it  follows  that  the  general  solution 
of  the  adjoint  system  can  be  v/ritten  as  a  linear  combination 
cf  the  13 J,  i.e. , 

(b)   U  =  Z  k •  u^     where  the  k.  are  arbitrary  constants. 

Substituting  (4)  into  (b)  results  in 

(6)  U  =  £k.  e"™^  (A 

Suppose  we  now  select  a  set  of  "special"  adjoint 

T 
splutions  U.  =  [un  .  uf, .  . ...  u  .  ]  ,         (j  =  l,2,...,n) 

which  meet  the  following  end  conditions: 

(7)  u   (T)   =  fi  ,   =/l  ,   i  =  j 

1J       1J    10  ,  i  /  3 

We  see  then  thst  the  terminal  values  x.(T)  of  the 
solutions  cf  the  original  system  cf  equations  can  be  deterg- 
ed by  successively  substituting  the  "special"  solutions 
of  the  adjoint  system  into  the  fundamental  formula  (3). 

That  is,  substituting  in  U.  in  (3),  we  get 

j 

_        T      rT 

a)     L  l  U..X    ]      =         c   Eu.^f.)  at 

iXJX^  Jo  i        J 

which,    since  u..(T)   =  1,   u..(T)   =  0,    i^j ,    reduces  to 

Jo  -'-J 

(9)        Xj(T)   =  I  ui(.(0)xi(0)   +f     (Lu.;.fi)dt,      (j=l,2,...,n). 


If  we  consider  U.  as  ;.  vector  in  a  cartesian  coordinate 

J 


set  with  base  vectors  ei?  then 
(10)   D,   =  I  Uije±. 


If  we  also  define  a  vector  X,  the  system  state  vector,  by 
the  equation 
(11)  X 


Ex.e. 


then  equation  (9)  can  be  written  as 

-  .    rT  -  - 

(12)  x.(T)   =   (U,.X)  +     (U..F)  at    (j  =  1,2,,.. ,n). 

c   Jo     J 


.] 


Since  the  IK  are  independent  solutions  of  the  adjoint 

J 

system  a  linear  combination  of  the  U.  is  also  a  solution, 
(12)  U  =  Ui  +  bAU2  +  ..,...+  bn_1Un 

where  the  b.  are  arbitrary  constants. 

A  solution  Li  defined  in  this  way  is  called  an  ad  .joint 
vector. 

From  (12)  we  have  the  following  equations,  obtained  by 
successively  substituting  in  the  U.  into  (3): 


rT     -   _ 

(  xA(T)   =   (u\.X)0  +     (UleF)  dt 

J   o 


x2(T)   =   (U2.X)   + 


-r  l     (U2.F)  dt 


(14)/ 


Un(T)   =   <VX)o  +     (VF)  dt 

J  o 
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If  we  form  any  linear  combination  of  the  above,  using 


say  the  b.  of  (13),  we  get 

J 

■T 

",F)dt 
o 


rT   -  ■ 

(15)   Xj.(T)  4  biXaCl)  +...+  bn  1xn(T)  =  (U.X)  +    (U.: 

n-i  n  o  J  0 


Equations  (14)  and  (15)  are  the  fundamental  equations 
for  the  scluticn  of  the  problem. 

We  shall  hereafter  consider  only  forcing  vectors  F  such 
that  the  scaler  predict  (U.F)  is  a  maximum.   In  the  next 
sections  it  will  be  seen  how  this  characterizes  the  solutions 
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4.   A  Theorem  on  Minimum  Time. 

In  this  section  it  will  be  shown  that  if  an  admissible 
trajectory  can  be  found  which  satisfies  a  maximizing  principle, 
then,  at  least  under  the  conditions  given  below,  it  furnishes 
the  desired  minimum  time,, 

An  admissible  trajectory  is  one  which  has  the  specified 
initial  and  terminal  points;  with  it  must  be  associated  an 
admissible  forcing  function .   An  admissible  forcing  function  F 
is  one  which  satisfies  the  constraints  on  its  components, 
|f . |  <  1.   In  this  section  no  specific  assumptions  are  made 
regarding  the  eigenvalues  of  the  matrix  A  in  the  system  of 
differential  equations 

(1)    X  =  A  X  +  F 
but  A  is  assumed  to  be  a  constant  matrix. 

a   point  of  static  stability  of  (1)  is  a  point  such  that 
an  admissible  force  can  be  chosen  to  constrain  the  system  to 

the  point.   That  is,  it  is  a  point  where   |  L  a. .x . |  <  1, 

j    J  ° 
(i  =  l,2,...,n),  and  x.  =  09  (i  =  l,2,...,n)  for  some  admis- 
sible F.   A  point  such  that  ;  Y.   a.-xJ  =1,  (i  =  l,2,...,n) 

j   J  J 
is  a  vertex  point  of  the  points  of  static  stability. 

Now  suppose  we  have  in  some  way  found  a  path  C  and  a 
vector  U  which  together  satisfy  the  following  relations: 

H.l.   The  path  C  is  admissible  for  some  value  of  T,  say  T  . 

H.2.   U  =  Uj  +  biUo  +  o00e+  b  -.IT   is  a  solution  to  the 
adjoint  system  of  differential  equations.   The  U  .  are  the 
"special"  solutions  defined  earlier,  and  the  b's  are  known 
constants. 
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H.3.   For  each  value  of  t,  the  force  F  associated  with  C 
is  such  that  the  scalar  product  (F  .U  )  is  a  maximum. 

H.4.  Each  component  u±*    (i  =  l,2,...,n)  of  U+  vanishes 
at  most  at  isolated  points  in  (0,T  ). 

Theorem.   The  trajectory  C  which  satisfies  H.l,  H.2,  H.3, 
and  H.4  yields  the  minimum  time  T+  and  is  unique,  compared 
with  all  paths  which  contain  a  point  of  static  stability 
that  is  not  a  vertex  point. 

Frocf .   Consider  equation  (3„15)  with  the  values  b1,...o,bn_1 

+ 
from  H.2  for  any  admissible  trajectory  for  time  T  . 

T+ 
(2)  Xlf  +  b,X2f  +....+  bn.Af   -    (U+.X)0  =j  Q    (if  .F+)dt 

We  shall  see  that  F+  maximizes  this  and  then  as  a  consequence 
T+  is  minimal.   Note  that  every  term  on  the  left  side  of  (2) 
is  determined  by  the  specified  conditions  X0,  Xf  and  the 
quantities  bi,...o,b  ^>  T  . 

Now  consider  the  force  F+  above  and  any  other  force  F1 


admissible  path  for  that  value  T  of  T. 

There  might  conceivably  be,  however,  a  path  C"  with  a 

smaller  value  T  of  T„   Assume  this  to  be  the  case  and  assume 

that  there  is  a  point  Xc  of  static  stability  on  C" . 

Let  T  -  T  =  AT,  and  let  t  be  the  time  associated  with  X„ 
e     '  s  s 

on  C". 

The  force  function  F'  constructed  as  follows  then  will 
lead  to  an  admissible  trajectory  C  with  T  =  T  : 

f."(t)  0  <  t  <  t 

J.  w> 

!■  la.x,       t  <  t  <  t  +AT 
"J  10  js         s       s 

fy(t-AT)  t  +AT  <  t  <  T+. 

The  middle  segment  of  C9  constrains  X  to  X  for  an 

interval  of  length  AT,  and  the  two  end  segments  then  lead 

from  X0  to  X  and  from  X  to  X^0 

Nov/  note  that 

T+  ts+AT 

(4)  f  [U+.(F+-  F»)J  dt   >  f     [U+.(F+-  F')]  dt, 

J  °  J  t 

s 

and,  for  at  least  one  value  of  i, 

tg+AT 

(5)  j  (f*  -  q)  u+  >  0 

ts 

since  ui  has  only  isolated  zeros  and  (f.-  f!)u.  =  (1  +  f!)|u. 

But  C  is  the  only  admissible  path  associated  with 
time  T  .   There  is  then  a  contradiction.   It  arises  from  the 
assumption  that  a  second  path  with  T  <  T+  existed.   This 
completes  the  proof;   if  such  a  path  is  found,  it  defines 
the  solution  and  is  unique. 
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Clearly  the  proof  holds  if  in  H03,  (F  ,U  )  is  a  minimum 
rather  than  a  maximum  on  cT  . 


Corollary.   If  the  matrix  A  has  real  distinct  eigenvalues, 
if  either  X0  or  X~  is  a  point  of  static  stability  which  is 
not  a  vertex  point,  and  if  a  curve  C  and  an  adjoint  vector 
have  been  found  which  satisfy  H.l,  H02,  H.3  with  none  of 
the  b's  of  H.2  zero,  then  C  defines  the  minimum  value  of  T 
and  is  unique . 

In  this  case  every  trajectory  contains  a  point  of  static 
stability  X0  or  X-  which  is  not  a  vertex  point.   Every  u. 
has  the  form  u.  =  Lc-.e~<3   with  at  least  one  c.  . 

distinct  from  zero  for  each  value  of  i.   It  has  at  most 

+ 
n-1  zeros,  and  U   satisfies  H04e   All  of  the  hypotheses  of 

the  theorem  are  satisfied  by  C  . 

If  the  values  of  ba , . . . . ,  b  -,,  T  were  known,  it  v/ould 

be  just  a  matter  of  integration  to  get  the  trajectory  by 

in- 
setting f .  =  sgn  u. 0   In  the  next  sections  these  will  be 

determined  by  a  method  of  successive  approximations. 

The  proof  of  the  theorem  in  this  section  is  due  to 

Professor  Faulkner . 
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5.   Example  of  a  Specific  becond-Order  Problem  and  Solution. 

In  this  section  the  solution  of  the  minimum- time  problem 
for  a  specific  second-order  stable  system  will  be  worked  out 
in  full  detail.   The  solution  consists  basically  of  the  fol- 
lowing steps: 

a.  The  adjoint  system  of  differential  equations  is  solved. 

b.  a   set  of  input  parameters  (X,b)  is  selected,  and  a  set 
of  "special"  adjoint  solutions  (with  the  corresponding  adjoint 
vector)  is  determined. 

Co   The  terminal  values  x(l),  y(T)  of  the  dynamic  system 
are  expressed  as  integrals  involving  the  "special"  adjoint 
solutions  and  the  vector  forcing  function,  and  a  trajectory 
is  determined  by  the  Maximum  Principle. 

d.   The  resulting  terminal  state  (x(T),y(T))  is  compared 
with  the  required  terminal  state  (Xf,Yf),  and  the  changes  in 
T  and  b  required  to  effect  (Xf,Yf)  are  determined  from  a  set 
of  differential-correction  equations. 

The  dynamic  system  which  will  be  considered  is  governed 
by  the  pair  of  linear  first-order  differential  equations 

9 

f  x  =   »  x+  y  +  f  |f  |  <  1 

(1)  <  . 

I  y  =  -2x  -  4y  +  g  |g|  <  1. 

The  system  of  differential  equations  adjoint  to  (1) 

can  be  written  directly  as 

• 

u  =  u  +  2v 
(2) 

v  =  -u  +  4v 
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i 


■ 


(3) 


The  characteristic  equation  of  matrix  A  is 
-1-m    1 
-2   -4-m 
The  eigenvalues  of  A,  roots  of  (3),  are  m  =  -2,-3. 


=  0 


Now,  if  mi,  m2  are  eigenvalues  of  matrix  A,  then  -m1,-m2 

T 
are  eigenvalues  of  -A  . 

2t   3t  T 

Hence  e   ,  e   and  the  eigenvalues  of  -A  will  define 

a  fundamental  set  of  solutions  of  the  adjoint  system  (2). 

Denote  these  solutions  by  U1,  U2  where 


(4)   U1   = 


U' 


V 


U; 


=   i  U' 


V 


Lv  j 


From  Appendix  I  we  know  that  for  each  distinct  eigen- 

T 
value  of  4  ,  the  corresponding  solution  of  the  adjoint 

system  can  be  written  in  the  form 


(5) 


U 


3   - 


■mit 


0L  C 


T      1 
where  -m.  is  an  eigenvalue  of  -A,  and  C^  is  the  corresponding 
J 

eigenvector. 

Since  the  eigenvalues  are  real  and  distinct,  we  know 
(see  Appendix  II)  that  the  components  of  CJ  are  in  the  same 
ratio  as  the  cofactors  of  the  elements  of  a  column  of  the 
characteristic  matrix  [A  -  m.l],   We  can  therefore  determine 
that 
(6) 


C1  =  [*cn 
c21 


=  12 


c2  =  r 


,c12   =   2  ! 

[c22j     2  J 

Then  the  solutions  U1,  U2  of  the  adjoint  system  are 
(7)   U1  =  e"mit  C1  =  e2tf2^ 

Uj 
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and 


3t  r«  1 


(8)  U2  =  e"m2t  C2  =  e~v 

L2j 

Since  U1,  U2  comprise  a  fundamental  set  for  the  adjoint 
system  (2),  IJ  =  kjU1  +  k2U2   is  also  a  solution,  i.e., 

(9)  U  =  kie: 


■2tr"2l     +     k2e3t   -2]       =[2     2irkie2tl 


2 


1      2 


l  .       3t 

i  k2e 


where  the  k's  are  arbitrary  constants, 

T 

We  will  denote  the  matrix  of  eigenvectors  of  -A  by  C, 

where 

(10)   C  =  [C1  C2]   =  f 2  2  " 

L1   2- 

-1 
Inverting  C  we  get  C   = ;  1  -1 


-  1 

~   2     J 


(ID 


Let  us  now  return  to   the  original   set  of  differential 

equations    (1)    which  we  may  write   as 

x+     x-     y  -  f     =     0 

y+2x+4y-g     =     0 

Suppose  we  now  multiply  the  first  equation  of  (11) 

by  u  and  the  second  by  v,  where  u^v  are  required  only  to 

be  continuous  functions  of  t  which  have  continuous  first 

derivatives.   Adding  the  two  products  and  integrating  the 

sum  from  0  to  T,  we  get  [4] 

T 
(12)     [u(x  +  x  -  y  -  f)  +  v(y  +  2x  +  4y  -  g)]dt  =  00 
J   o 

Note  that  if  x,y  are  solutions  of  the  original  set  of  equations, 

then  u,v  used  in  this  way  are  called  Lagrange  multipliers. 
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Integrating  (12)  by  parts  we  get 

(13)        [ux  +  vy]      +/        [x(-u  +  u  +  2v)    +  y(-v   -  u  +  4v)]dt 


J   o 

■/' 


[uf  +  vg]  dt  =  0. 


If  we  now  choose  u?v  as  solutions  of  the  adjoint  system, 
we  have 

f  -u  +  u  +  2v  =0 


(14) 


-v  -  u  +  4v  =  0 


If  u,v  are  thus  chosen,  (13)  reduces  to 
T     r  T 

o 


T     r  T 
(15)    [ux  +  vy]o   =  /    [uf  +  vg]. 

-J    o 


Note  that  (15)  is  a  particular  example  of  the  Fundamental 
Adjoint  Formula  of  section  20 

Suppose  we  now  choose  a  set  of  "special"  solutions  of 
the  adjoint  system  as  follows s 


(i6)    Ui  =  rur 


vi 


and 

(17)   U2 


^2 


L 


such  that  Ui(T)  =  1 

v±(T)  =  0 

such  that  u2(T)  =  0 

v2(T)  =  1 


Substituting  Uj.  into  (9)  we  get,  evaluating  at  t  =  T 


(18)   Ui(T)  =  |Ul(T) 

va(T) 


kne^1]  H'l] 


k21e3TJ    0 


=  j  2  2 

L1  2, 
Note;   The  second  subscripts  of  k11}  k21  indicate  that  they 

correspond  to  the  first  solution  Ulo 
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Multiplying    (18)    through  by  C"     we   get 


(19) 


r.  2T")      _      P-l      , 

i  kue  C 


k21e 


3T 


L°J 


where  the  cj  ,  are  elements  of  C 

-2T 

Therefore  kn  =  e 


r « 5  ' 

j  Cll 
;  C2i 

=1 


-21 


=  e 


■3T 


>2 


1 
=1/2 


Then 

(20)   Ui(t)] 

[    Vi(t) 


=  ;  2    2 
1    2 


2(t-T) 


L  =-2 


3(t-T) 


Similarly,  substituting  U2  into  (9)?  we  determine  that 
(21)   k12  =   -  e"2T 

k22  = 
so  that 


■3T 


(22)   u2(t)  I  =  [2     2 


r    2(t-T)i 


v2(t)  i    1  2  ;  e3(t°Tj; 

Substituting  Uj,  into  (15)}we  get 

T 
(23)   x(T)  =  u1(0)x(0)  +  V!(0)y(0)  +  I   [u^  +  vlg]dt. 

J   o 


Substituting  U2  into  (15),  we  get 


,T 


(24)  y(T)  =  u2(0)x(0)  +  v2(0)y(0)  +1   [u2f  +  v2g]dt. 

J  o 

Let  us  now  consider  the  "special"  adjoint  solutions 
as  vector  functions 

(25)  Uj.  =  Uiij  +  Vjez 
U2  =  u2ei  +  v2e2 

where  el9   e"2  are  the  base  vectors  for  a  cartesian  coordinate 

set  in  two  dimensions, 
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Since  Uj  and  U2  are  independent  solutions  of  the  adjoint 
system,  a  linear  combination  of  these  is  also  a  solution,  i.e., 

(27)  U  =  Ui  +  bU2 

=   (Ui  +  bu2)e!  +  (Vi  +  bv2)e2 

=  u  ex     +  v  e2 
Combining  (23)  and  (24)  in  the  same  linear  combination,  we 
obtain 

(28)  x(T)  +  b  y(T)  =  (u!+bu2 )Qx(0)  +  (v1+bv2)oy(0) 

T 
+    [(u1+bu2)f  +  (v!+bv2)g]dt. 
J  o 

which  can  be  written  as 

T 

(29)  x(T)  +  b  y(T)   =  u(0)x(0)  +  v(0)y(0)  +    [uf  +  vg]dt 

■   o 

Defining  vector  R  =  xei   +  ye2,   we   can  rewrite (29)    as 

T 

(30)  x(T)   +  b  y(T)      =      (U.R)        +  (U.F)dt. 

°  o 

where  F  =  f  ej.     +     ge2. 

Suppose  now  that  we  desire  to  drive  the  system  from 
(x(0),y(0))  to  (Xf?Yf)  in  minimum  time.   As  we  have  seen 
earlier  (section  4),  this  may  be  accomplished  by  applying 
Pontryagin's  Maximum  Principle  to  (30),  i.e.,  by  choosing 
the  vector  F  =  F  such  that  the  projection  of  F  on  the 
adjoint  vector  U  is  at  all  times  a  maximum.   That  is, 
we  choose  F  such  that  its  components  satisfy  the  following 
relations? 

ft       =   sgn  u 

(31)  /   + 

I  g   =  sgn  v0 
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Equation  (30)  then  can  be  rewritten  as 

T  .  _+ 

(32)  x(T)  +  b  y(T)   =   (U.R)n  +      (U.F  )dt 

0     o 

while  (23)  can  be  rewritten  as 

_  _       T  „  „+ 

(33)  x(T)   =   (lil0H)o  +    (Ui.F  )dt  , 

-  o 

and  (24)  can  be  rewritten  as 

-  -      ,  ?  _  -+ 

(34)  y(T)   =   (U2eR)0  +    (U2eF  )dt. 

o 

However,  the  value  of  T  which  we  have  selected  is  only 
a  "guess"  at  the  actual  minimum  time  required  to  drive  the 
system  to  (Xf,Yf),  and  the  constant  b  is  only  a  "guess"  for 
the  linear  combination  of  U1?U2  to  determine  the  proper 
adjoint  vector  associated  with  the  desired  terminal  state. 
We  must  therefore  compare  the  resulting  state  (x(T),y(T)) 
with  the  desired  terminal  state  (X~,Y~)  to  determine  the 
errors 

(35)  J  Xf  -  x(T)   =  AX 
]  Yf  -  y(T)   =  AY 

and  determine  the  changes  $T$     #b  required  in  the  input 
parameters  to  drive  the  system  to  the  desired  terminal  state 
For  this  we  will  derive  a  set  of  equations  to  determine 
Ol,  6b0   Suppose  we  replace  T  and  b  by  nearby  values 

(36)  f  T ^  T  +  6T 

b  — — *►  b  +  6b 


Now,  writing  x(T)  =  x(T9b)  and  y(T)  =  y(T,b)   since  the 

terminal  values  are  implicit  functions  of  b,  we  get 

(37)  (   x(T,b)  *>  x(T+  ST,b+  6b)  =  x(T,b)  +  &X 

I  y(T,b) *-  y(T+  6T?b+  db)  =  y(T,b)  +  6Y 

plus  terms  of  higher  order „ 
Now 


(38) 


r   £*  -    5X 

6T 
6Y  =   5Y 


&T 


&T 


6X 


Sb 


db 
^6b 


dT        5b 

where  <SX,  6Y  denote  the  total  differentials  of  x(T,b) 
and  y(T,b)  respectively,  i0e0? 


(39)  < 


-£  x(T)  =  & 

dT  K    J  6T 

-*  y(T)  =  ^ 

6T  yK    J  £T 


_P_     /rnv  5X 

5b  XU;  5b 

_5   m  5Y 

db  yki;  6b 


We  can  rewrite  (38)  as  a  single  matrix  equation 


(40)   [*6x'i 


i  ax  ex 

|  ST  6b 

5Y  5Y 

.61  6b 


'*T 
6b 


6b 


where  £"  is  the  matrix  defined  in  (40) „   This  matrix  is 
called  the  "correction"  matrix0 

Multiplying  (40)  through  by  £"  "   and  setting  &X  =  AX 
and  6Y  =  AY,  we  get 


(41) 


ST 

6b 


=  € 


•1 


6x1 

6Y 


Xf-x(T)j 
Yf.»y(T)| 


Equation  (41)  will  be  referred  to  as  the  "correction  equation. 
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Hence  we  have  determined  the  changes  ST,  6b  in  the 
input  parameters  required  to  drive  the  system  to  (X~ ,Yf ) ,  so 
that 

(42)  f    x(T  +  £T,  b  +  6b)   =  Xf 
1  y(T  +  6T,  b  +6b)   =  Yf 

except  for  higher-order  terms 0 

Therefore  we  have  solved  the  minimum-time  problem 
assuming  we  can  determine  the  elements  of  the  "correction" 
matrix,  and  if  the  routine  converges „ 

We  will  first  determine  dX  . 

aT 
-  -     rT  -  -+ 

where  X  =  x(T)  =   (U^R)   +     (U^F  )dt. 

°     Jo 

On  a  minimum-time  trajectory  each  component  of  a 

two-dimensional  forcing  vector  will  switch  (change  sign) 

at  most  one  time  if  the  eigenvalues  are  real.   These 

switching  times  are  defined  by  the  equations 

(   u(tj)   =  0  9    f  switches  from  +1  to  +1  at  t  =  tle 

(43)  <^ 

I   v(t2)      =     0    ?  g    switches   from  +1  to  +1   at   t   =  t2. 

Order   tlf   t2      as   tk,  0  <  k  <  K  <  2 

Define  (  0     =  t0  k  =  0 

1  T     =   tK+1  k  =  K 

so   that     0  <  k  <  K  <  2. 

rtith  the   above  notation,   we   can  express    (33)    as 

K  t 

(44)  x(T)      =      (D 


r1.B)0  +  i    fk+1   GiTHt 

k=0  J  t,. 


'k 

where  the  integrand  is  continuous  over  each  subinterval 

cf  integration  t,  to  t,  m. 
fc        k     k+1 

24 


To  simplify  notation  we  will  make  the  following 
definition: 

(45)  Ik  =   (UlF*)  ,         tk  <  t  <  tk+1 

Then  (44)  can  be  rewritten  as 

K  fVt-l 

(46)  X  =   (U1oR)0  +  Y.   J  Ik  dt- 

o   tk 

Applying  Leibniz'  Rule  for  the  differentiation  of  a  definite 
integral  the  limits  of  which  are  functions  of  a  parameter  [5J, 
we  obtain 

(47)   3  =  4x^o +  y.k+1  fik  dt   +  i  vw  ew 

-    I    ik(t  )  etk 

or  if  [   ].   denotes  the  discontinuity  (jump)  in  a  function 
^k 

at  point  tk,  we  can  write  (47)  as 

(48)  ||  „  &.£)      ♦  fT  (fSi.Ddt  ♦  (L-..^), 

J  o 

-  E  CtJi.LF*"]).   6tk 

i       \  eT 

where  K  =  0,  1,  or  20 

Now,  for  fixed  b,  it  is  readily  seen  that  £tk  _  1 

q1     -  l,  so 

so  that  (48)  reduces  to 

(49)  |X     =      (gi.S)o  +^"    cffli.Ddt      +(U1.F+)T 

rk 
1 
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Similarly, 

(50)   ff  =  ^o  +f  ^».F+)dt  ♦  a2.F+), 


-  E  (U2.[F+])t 

1         ^k 

By  the  same  methods  we  can  determine  that  eX  and  eY 

6b      eb 


can  be  expressed  as 

K-l 

(bi)   ex 

eb 


k=0    1 


k+l 


et 
db 


k+l 


and 


K-l 


(52)   jX  s  _  £   (U.[F+]) 


eb 


k=0 


t,  ^    dtk+l 

k+1  eb" 


ft  Vv  ft  Y 

To  completely  determine  rr  end  rr- 


ft  t 

we  must  determine  , k+l 

eb 


These  can  be  determined  by  making  the  following  linear 
approximation  of  the  change  in  the  switching  times,  denoted 
by  6ti   (and  6t2)s 

eu 

tl 


(53)    fit!   = 


on 


L(gw  +  |M,b)]ti 


u 


and  since  u  =  u(T-t),  it  follows  that  u  = 


9u 
6T 


so  that  (53)  reduces  to 
(54)    6tA  =  6t 


eu 
e_b. 


6b  =  6t 


6b 


so  that 

(55)    at!  _ 

eb 

Similarly, 

u 

• 

(56)    eta  - 

eb 

V 

• 
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Having  determined  6t^+-,   we  can  now  write  the  final 

eb 


expressions  of  e_A  and  e_Y  as 
eb     eb 


(57)   ex 

5b 


UiLf+]  -u2 


u 


+  I vi[g  J  j-v2 


v  /J  t. 


and 


(58)   eY  _ 

eb 


u2[f  ]  /  -u2 


U 


ti     L      \  v  /Jt2J 

We  have  now  determined  all  the  elements  of  the  "correction" 
matrix,  so  that  the  minimum-time  problem  is  solved,  if  the 
routine  converge s0 

In  this  section  we  have  developed  the  solution  of  the 
minimum-time  problem  for  a  specific  second-order  system„   However 
each  step  of  the  solution  has  been  derived  in  general  notation 
before  substituting  in  the  numerical  values.   It  is  therefore 
unnecessary  to  repeat  the  development  for  the  general  case  for 
second-order  systems.   In  the  succeeding  sections,  the  solution 
of  the  problem  for  the  general  case  will  be  developed  for 
higher-order  systems. 

An  interesting  property  of  second -order  systems  was 
discovered;  it  is  the  f ollcwing i   If  a  minimum-time  trajectory 
for  second-order  systems  has  two  switching  times,  the  time 
between  these  switching  times  is  an  invariant  of  the  dynamic 
system.   An  elementary  proof  of  this  result  is  given  in 
section  8. 
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6.   Derivation  of  Solution  for  Third-Order  System. 

Suppose  we  have  a  set  of  three  inhomogenecus  linear 
first-order  differential  equations  with  constant  coefficients 
representing  a  dynamic  system 


(1)   X  =  A  X  +  F 


iT 


where  X  =  [xa  x2  x3]x9   F  =  [fa  f2  f3]x,   and  A  is  a  3  x  3 
matrix  with  constant  elements  a.,  „   The  eigenvalues  of  A 
are  assumed  to  be  real  and  distinct 0 

The  characteristic  equation  associated  with  matrix  A 
is 

(2)  |  A  -  ml  |   =  0, 

1  •  c  •  • 

(3)  an-m   B-i2  ai3 

a21   a22~m    a23 
a31     a32     a33"m 

which,  when  expanded  by  determinants,  becomes  a  cubic 


equation  in  m  of  the  form 
(4)   m  +  ba.ni  +  b2m  +  b< 


0. 


Denote  by  ml9   m2 ,  m3  the  roots  of  (4),  i.e.,  the  eigenvalues 

of  matrix  A0 

T 
Then  -mj ,  -m2 ,  ~m3  are  eigenvalues  of  4  . 

Therefore,  if  we  consider  the  system  adjoint  to  (1), 

which  is 


(5)   U  =  ~AX   U 


where  U  =  [ux   u2  u3] 


T 


then  e 


■mat   -m2t   -m3t 


T 


and  the  eigenvectors  of  -A  will 


define  a  fundamental  set  of  solutions  of  the  adjoint  system, 
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Let  us  consider  the  solution  of  the  adjoint  system  (5) 
in  terms  of  the  eigenvectors  of  the  matrix  of  coefficients, -A 
Ae   have 

(5)   u  = 


T 


T 


Assume  a  solution  of  the  form 

(6)  U  =  e   C1   where  C'^  [cj.  c2  c3]   is  a  constant 
vector c 

Differentiating  (6)  with  respect  to  t,  we  obtain 

(7)  u  =  r  ert  c1 


T 


Substitution  of  (7)  and  (6)  back  into  (5)  gives 

(8)  re   c  =  -A  e   c? 

rt 
and  since  e   is  a  scalar  quantity  appearing  on  both  sides 

of  the  equation,  it  follows  that 

(9)  r  C1  =  -AT  c\ 
Therefore,  if  we  choose  r  as  an  eigenvalue,  -m.,  of  -A" 

then  C'is  by  definition  the  corresponding  eigenvector,,   Hence, 

T 
for  each  eigenvalue  of  -A  we  have  a  unique  solution  of  the 

adjoint  system  if  the  eigenvalues  are  distinct,,   Therefore, 

we  have  the  following  fundamental  set  of  solutions  of  (5) s 

u2 


-*-  c 

■mot 


2^  Q2 


0 


3 


■m« 


Then  any  linear  combination  of  U1,  U2,  U3  is  also  a  solution 

(10)   U  =  kie^mit  C1  +  k2e=m2t  C2  +  k3e"m3t  C3 
where  the  k's  are  arbitrary  constants.  If  (10)  is  written 


in  expanded  form,  we  get 


(11) 


T7      i    -mit  r 

Ij  =  kje  x        en 

C21 


koe 


=mPt 


j  ^21 
[C31 


Cl2 
C22 
C32 


koe 


=m»t 


I  c13 
c23 
C33 
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Equation  (11)  can  then  be  expressed  as  the  product  of 
two  matrices,  as  follows: 


Hi 

™ 

u3 

en   c12   c13 
Cj31  c22  C23 

C31   C32   C33 


m  1 1  -, 


1-1 
-nigt 

-m3t 


1 
k2e 
JL  k3e 


where 


C  = 


C11   c12   c13 

C21    c22    c23 

c31   c32   c33 


is  the  matrix  of  eigenvectors  of  =A  . 

Then,  the  general  solution  of  the  adjoint  system 
can  be  written  as 


(12) 


11 ! 

^2 

u3 

Hi                             -J 

,   -nut  -, 
kie   x 


k-e 


-mot 


Lkae-"1"*. 


We  will  now  chocse  a  "specialnsolution  Ui  =  [uu  TJ2i  u3i] 
such  that  Un(T)  =  1,  u21(T)  =  09  u31(T)  =  0„   Substitution 
of  Ui  into  (12)  evaluated  at  t  =  T  gives  us 


(13) 


Un(T) 
u2i(T) 
Usi(T) 


=  C 


kne 

k2ie 

k3ie 


•mil  - 
-m2T 

■m3T 


where  the  second  subscripts  of  the  k!s  indicate  constants 
corresponding  to  solution  Ul9 

Multiplying  (13)  through  by  C"  we  get 


(14) 


r,    -nuT" 
I  kne  x 

-m2T 


k2ie 
k3ie 


=m,T 


1 

1 

Oil 

0 

= 

1 

C21 

0_ 

c3i 

where  the  c! .   are  elements  of  C' 
J-  j 


1.   See  Appendix  II  for  determination  of  components  of  C^ 

in  terms  of  cofactcrs  of  [a  =  m»l]0 

J 


$0 


Then 


(15) 


•    m«T 
=  en  e  1 

1    nioT 
=  c21  e  * 

«    m3T 
=  c3i  e  3 


Substitution  of  (15)  back  into  (12)  gives 

■n^t 


(16) 


"11 
^31 


•11' 


k21e 


=m2t 


-81' 


=  c 


cl,  em*(*-t) 
c21  e**^ 


Similarly,  we  can  choose  a  second  "special"  solution 

T 
U2  =  Lui2  U22  ^32]   °f  "tne  adjoint  system  satisfying  the 

end  conditions  u12(T)  =  0,  u22(T)  =  1,  u32(T)  =  0,  which 

results  in 


(17) 


Ui2 

U22 

Ju32_ 

and,  choosing  a  third  "special"  adjoint  solution 

U3  =  [U13  u23  u.33]   satisfying  the  end  conditions 

u13(T)  =  0,  u23(T)  =  0,  u33(T)  =  1.  we  obtain 

■  •   iMT-t)] 


(If) 


U23 
^33 


c13  e 

1   m2(T-t) 

c23  e 


If  we  group  (16) 
equation,  we  obtain 


m3(T-t) 
_c33  e 

(17) ,  and  (IS)  into  a  single  matrix 


(19) 


ull  u12  ^13 

u21  u22  u23 
u31  ^32  u33 


=  c 


(20) 


(21) 


[U]  =  c 


[U] 


c'21e^(T-t) 

>i(T~t) 


c28em^T"^ 
c*2em*(T-t) 


c'13erai(T"^ 
c.9em2(T-t) 

c,3em3(T»t) 


0 
0 
=1 


0 
,m2(T-t) 

0 


0 
m8(T-t) 


—  C  E  C 

where  E  is  the  diagonal  matrix  with  jthdiagonal  element  emJ^T"t^. 
and  [U]  =  [Uj  U2  U3]. 
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Since  Ui,  u2,  13 s  are  solutions  of  the  adjoint  system 
then  &ny  linear  combination  of  these  is  also  a  solution: 

(22)   U  =  Ui  +  bjUg  +  b2U3 

T 

where  U  =  [uj  u2  u3]  j 

i0e0,  [  ua  =  Un  +  biU12  +  b2ui3 
I   u2     =     u21  +  biu22  +  b2u23 
u3  =  u31  +  b!U32  +  b2u33. 
Written  in  matrix  form?  the  above  set  of  equations  becomes 


(23) 


Ui       I  Un  u12  u13 

u2   =  I  u21  u22  u23 
u3       I  u31  U32  U33 


1 

bj 

b2 


Let  e1?  e2,  e3  represent  the  base  vectors  of  a  three- 
dimensional  cartesian  coordinate  system*  We  can  then 
express  Ui,  U2?  b3  as  the  vectors 

£1   =  Une!  +  u21e2  +  u31e3 

(24)  U2  =  u12e!  +  u22e2  +  u32e3 
U3  =  u13ei  +  u23e2  +  u33e3 

Then  U  =  UA  +  b1U2  +  b2U3  can  be  written  as 

(25)  U  =  (uu+  bxu12+  b2u13)e1  +  (u21+  biU22  +  b2u23)e2 

+  (u31+  biU32  +  b2u33)e3 
which  can  also  be  written  as 

(26)  U  =  i^ei  +  u2e2  +  u3e3 

A  solution  U  obtained  in  this  way  is  called  an  ad  .joint 
vector o 

If  in  the  Fundamental  Adjoint  Formula  we  let  t±    -   0 
and  t2  =  T,  we  have,  from  (2.11), 

T 


(27)   [  I  uiXi]o   =  J         (  Zu  f  )dt9 
1  o     i 
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If  we  desire  to  determine  x.(T),  we  chocse  the  adjoint 
solution  U  =  U..   This  reduces  (27)  to 

(28)  x.(T)   =  (L  uijx.)o  +J   (Iu..f.)dt   (3=1,2,3) 
which  can  be  written  in  equivalent  vector  notation  as 

-   -  rT   -   - 

(29)  X;J(T)   =   (Uj.X)Q  +  J    (Uj.F)dt,         (3=1,2,3), 

where  X  =  x^j  +  x2e2  +  x3e3° 


we  wish  to  drive  the  system  from  state  X(0)  to  X~  in 
minimum  time  T.   We  therefore  consider  a  linear  combination 
of  the  x.(l),  using  say  the  constants  b.  of  (22).   We  get 


(30)    Xl(T)  +  biX2(T)  +  b2x3(T)  =  [ (Uj  +  bxU2  +  b2U3)„X]0 

rI     _  _  _       _ 

+  /     Ltti  +  biU8  +  b2U3).F]dt 
J   o 


which  can  be  written  as 

rT     _. 

(31)  Xi(T)   +  biX2(T)   +  b2x3(T)    =      (13.X)Q  +  (U.F)dt. 

J   c 

Applying  Pontryagin's  Maximum  Principle,  i.e.,  by  choosing 

F  =  F  such  that  f .  =  sgn  u.,  we  obtain 

T 

(32)  Xi(T)  +  blX2(T)  +  b2x3(T)  =   (U.X)Q  +  f   (U.F+)dt 

We  now  compare  the  resulting  terminal  state  X(T)  with 
the  desired  terminal  state  X~,  and  define  the  component  errors 
as 

'   Xlf  -  Xl(T)   =  AX! 

(33)  <  X2f  -  x2(T)   =  AX2 

^3f  "  Xa      =  ^*-3  i 
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If  any  AX.  (j  -  1,2,3)  is  not  zero,  we  must  determine 
the  changes  <5T,  <5bi,  6b2   in  the  input  parameters  required 
tc  effect  X*..   Suppose  the  initial  input  parameters  are 
replaced  by  nearby  values 

'  i  >  t  +  6T 

bj  3-  bA+  6bj 


(34) 


bo  - 


■*-  bp+  6bp  „ 


Then  the  resulting  changes  6X.  in  the  terminal  values  x.(T) 
are  determined  to  be  (for  j  =  1,2,3) 

T 


Ka 


\Mi.  x)_  +  f     (j|Ej.F+)dt-I     ._  (u.-Lf^]).1 
I  ei         o    j  o    si  i=i  k=1    ij     i    tu 


lit 


v1  l^^l 


i=l  k=l 


u. 


'ik 


+   f*(T)  ]     6T 


6b! 


_____ — _. 


ui 


ik 


6b. 


where  K.  denotes  the  number  of  zeros  of  u. »  and  t..  the  cor- 

+  1         irl 

responding  times.   Equation  (35)  has  the  form 

(36)  6XA      =  f&j  6T     |$j  6b,   +  5|J  Ob2     (3=1,2,3) 


6T 


5bi 


ab: 


Setting  oX.  =  -X. ,  and  inverting  the  "correction" 

j     j 

matrix  defined  by  (36),  we  can  determine  the  required  changes 
6T,  Obj  ,  6b2   to  effect  Xf0 


i.   (^[fpu 


ik 


=  u,  ,fT 

ij  i 


tik+o 


1J*- 


t.,  -0 

ik 
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Derivation  of  Solution  for  I\l-th  Order   Systems. 


th 


Let  us  consider  a  linear  dynamic  system  of  the  n   order 
represented  by  the  system  of  differential  equations 

(1)  X  =  A  X  +  F 

where  A  is  an  n  x  n  matrix  with  constant  elements  a.,  and 

T 

with  real  and  distinct  eigenvalues.   X  =  [xA  x2  ....  xl 

and  F  =  [fj.  f2  °  = .  *  fn]T» 

The  system  adjoint  to  (1)  can  be  written  directly  as 

(2)  U  =  -  aT  U 

where  U  =  [uj  u2  <> .  .  3  u  ]  .   We  know  from  Appendix  I  that 


equation  (2)  has  solutions  of  the  form 
(3) 


U  =  Ek4  e'm^   C^ 


T  i 

where  -m.  are  eigenvalues  of  matrix  -A  ,  and  the  C°  are  the 

corresponding  eigenvectors. 

Equation  (3)  can  also  be  written  in  matrix  form 

as 


(4) 


U  =  C 


'l^e'*^ 


ko.e 


-mot 


^.-n* 


where  C  is  the  n  x  n  matrix  composed  of  the  eigenvectors  CJ , 
i.e.,  C  =  [C1  C2  .  .  „  0  C  ]  where  the  eigenvectors  CJ  are  n  x  1 

column  matrices o 

iT 
If  v/e  choose  a  "special1"  solution  Ui  =  [u^  u21  ,,.u  J 

to  the  adjoint  such  that  un(T)  =  1,  u21(T)  =  . ...  =  u  -,(T)=0, 

then  equation  (4)  evaluated  at  t  =  T  becomes 
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(5) 


'un(T) 

u81(T) 


unl(T) 


=  C 


klle-^" 

1 

k21e"^T 

= 

0 

Knle 

0 

_     - 

1 


Multiplying  (5)  through  by  C   we  get 


(6) 


i    -mil 
kne   1 

k21e   2 


knle 


•mnT 


[1 

•1  0 


L°J 


1 
Cll 

C21 

_°nl. 

where  the  c.  .   are  elements  cf  C 

Therefore,  the  constants  k.-,  corresponding  to  solution 

U!  are  determined  to  be  as  follows: 

.mjT 


k 


1 1 


(7) 


kc 


=  c21e 


•  _rn2T 


k 


nl 


1   innT 
c  t  e  n 
nl 


so  that,  substituting  (7)  back  into  (4)  we  get 


(8) 


Un(t) 

u21(t) 


unl(t) 


=  C 


kiie-"1 


k2ie 


nijt' 
m2t 


k  e°mnt 
Knle 


=  C 


c.iem2(T-t) 


.   mn(T-t) 
cnle 


so 


Therefore,  if  in  general  we  choose  the  "special" 

T 

lutions  U.  =  [u1-  u2jeo* tUn1^   to  the  adJ°int  syst 


em 


such  that  u. .(T)  =  & 


l  0   i 


=  J 
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we  get 


(9) 


u^Ct) 
u2j(t) 


V*> 


=  c 


,   m2(T-t) 


'   mn(T-t) 


so  that,  by  grouping  the  n  sets  of  solutions  so  obtained 
into  a  single  matrix  equation,  we  obtain 

>!<*-*>    0     .  .  . 


(10) 


un  u12 

^21  ^22 


unl  un2 


"in 
u2n 


u. 


nn 


=  C 


0 


0 


m2(T-t) 


0 


0 
0 


mn(T-t) 


•1 


Since  the  U»  are  independent  solutions  to  the  adjoint 
j 


system,  a  linear  combination  11=11!+  biUg  + 


u„, 


^n-l"n: 

where  the  b's  are  arbitrary  constants,  is  also  a  solution. 
A  solution  obtained  in  this  manner  is  called  an  adjoint 
vector. 

If  e1?  e2?»°°°en  represent  the  base  vectors  for  an 
n-dimensiorial  cartesian  coordinate  set,  we  can  express  the 
adjoint  vector  U  as 
(11)   5  =  Y.    u±   e± 


where 

(12)        u±  -     uil  +  b!Ui2  + 


©      O      Q      © 


.       +     b  nll. 

n-lin 


\  1       "*"       I  joi    s    o    o    qII/    e 


For  a  minimum- time  path  we  choose  F  =  F  such  that  the 


components  of  F  satisfy  the  relation 
(13)   fT  =   sgn  u, 

1  1 


\j.    i«£j«  o  «  «  «nj  o 
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From  the  Fundamental  Adjoint  Formula  v/e  have 

T  r  T 

(14)  [  £  u.x.]        =  (I  u   f+)    dt 

iiic  Jo         ill 

which,  upon  successive  substitution  of  the  "special"  adjoint 

solutions  U.  reduces  to 
j 

rT  _    , 

(15)  x,(T)   =   Lu   (0)x.(0)   +     (  zu  f  )  dt 

J         i^  ^  o   i 

for  j  =  l,2,o.,n„   Equation  (15)  can  alsc  be  written  in 
vector  notation  as 

(16)  x.j(T)   =   (Uj.X)Q  +/   (Uj.F  )  dt    (j  =  1,2,., ,n). 

Again,  as  in  the  previous  sections,  we  desire  to  drive 
the  system  from  X(0)  to  X~  in  minimum  time  T.   We  therefore 
consider  a  linear  combination  of  the  x.(T),  using  the  constants 
previously  described  in  defining  the  adjoint  vector.   We  have 

T . 

(17)  xi(T)  +  biXgCT)  +  ...+  bnalxn(T)   =  (U.X)   +    (U.F  )dt 

J    o 

where  U  =  lij.  +  bjL2  +  . ..  +  kn_i^n  is  the  adjoint  vector. 

Comparing  the  resulting  terminal  state  X('i')  with  the 
desired  terminal  state  Xf ,  we  define  the  component  errors 
as 

(18)  AX..   =  Xjf  -  Xj(T)  (j  =  1,2,..., n). 

If  any  AX.  is  not  zero  we  must  determine  the  changes  in  the 
input  parameters  required  to  effect  X-. 
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Suppose  the  initial  values  of  the  input  parameters  are 
replaced  by  nearby  values 
T »  T  +  5T 


(19) 


bi  +  &bx 


Vl— *Vl+  6bn-r 
Then  the  changes  <$X.  resulting  in  the  terminal  values  x.(T) 
can  be  expressed  as  (fcr  j  =  l,2,.„.,n) 

(20)   6X.  =   (§2o.X)c  ^j*1  (§§J.F+)dt  +  f+(T) 


i=l  k=l 


n   La 

i=l  k=l     1J    x   t 


cT 


I   I1    (V!!!!!!1 


u- 


6bx 


'±kJ 


ik 


i=l  k=l 


u- 


'ik- 


6  b. 


n 


„"i 


Equation  (20)  has  the  form  of  the  total  differential 

6X 


(21)   <5x 


tf3*T 


ebi   x 


•  •  •  +  ft3  *  V 

n-1 


for  (j=l,2,...,n). 

Setting  AX.   =  OX.  and  inverting  the  "correction" 

J  J 

matrix  defined  by  (21) ,  we  can  determine  the  changes 

Si',  6bi,...„,  6b  n  required  to  effect  X~,  if  the  routine 

converges. 
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£„   Conclusions  and  Comments, 

The  solution  of  the  minimum-time  problem  for  second  and 
third -order  systems  was  programmed  and  run  on  the  CDC-1604 
computer,  using  the  Runge-kutta-Gill  method  of  numerical 
integration  to  calculate  the  required  integrals.   It  has  a 
desirable  feature  that  it  is  easy  to  vary  the  time  intervals 
any  time  during  a  run„ 

The  convergence  of  the  differential-correction  routine 
appears  to  be  fairly  critical,  especially  for  stable  systems. 
That  is,  if  for  a  given  set  of  input  parameters  the  error 
X~  -  X(T)  is  toe  large,  the  routine  will  not  converge. 

The  switching  times  ti?  t2  for  second -order  systems  were 
determined  exactly,  so  that  the  degree  of  convergence  was 
limited  only  by  the  degree  of  precision  provided  in  the  numerical 
integration  (limited  by  the  size  of  At)  and  other  routines. 

Convergence  of  the  routine  was  obtained  to  an  accuracy  of  better 

-12 

than  six  decimal  places  (AR2  =10   ). 

For  the  third-order  systems  the  switching  times  were  not 
determined  exactly  §   if  a  component  of  the  optimum  forcing 
function  changed  sign  in  a  time  interval  At,  it  was  assumed  to 
have  changed  sign  at  the  end-point  of  the  interval.   In  this 
case,  convergence  of  the  routine  was  limited  to  about  three 
decimal  places  (AR2  -  10   )  for  At  =  .005.   Convergence  to 
greater  accuracy  could  be  achieved  by  reducing  the  size  cf  At 
in  the  numerical  integration  routine. 

The  method  of  solution  appears  to  be  applicable  to  more 

general  problems „      The  case  where  the  end -point  Xf  varies  with 

time  T  may  be  handled  by  incorporating  into  the  correction 
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equations  (701£)  a  term  invclving  dXj_/d.T,  as  was  dene  by 
Faulkner  and  wvard  [10]  in  a  rocket  trajectory  problem. 

The  solution  curves  have  several  interesting,  properties. 
As  noted  in  section  5,  if  the  solution  for  a  second -order 
system  has  two  switching  times,  the  difference  is  a  constant 
of  the  system0   A  proof  is  given  here  that  this  is  generally 
the  cases 

Let  us  consider  a  linear  dynamic  system  represented  by 
the  system  of  differential  equations 

(1)   X  =  A  X  +  F, 
which  can  be  written  in  expanded  form  as 

x  =  ax  +  by  +  f  |f |  <  1 

y  =  ex  +  dy  +  g  |g|  <  1 

where  the  elements  a?b,c,d  are  constants,  and  the  eigenvalues 
of  A  are  real  and  distinct „ 

The  eigenvalues  of  matrix  A  (rccts  of  the  characteristic 
equation  |A  -  ml |  =0  )  are 


(2) 


(3) 

where 


m- 


a  +  _d_  +  r 
2 


m2 


a  •+  d  -  r 
2 


(4)   r  =  /(a=d)2  +  4b  c   „ 

The  "special"  solutions  Ul3  Ue  (described  in  earlier 
sections)  of  the  system  adjoint  to  (1)  can  be  expressed  jointly 
in  the  following  matrix  equation: 


(5)    [13!  U2] 


xxt      u2 

V^    Vo 


=   C 


em!(T-t)    Q 


m2(T-t) 


where  C  is  the  matrix  of  eigenvectors  of  -A  , 
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The  general  solution  h   =  Uj  +  kU2  of  the  adjoint  system 
can  then  be  expressed  as 


(6)        U      =   ju  Uj      u2 

vi      v2 


i: 


l' 

k 


l0r9*iW~t) 


\     L 


o 


0 

m2(T-t) 


c 


rl"> 

k 


Let  t1?t2  denote  the  switching  times  of  f  and  g  respectively, 

i.e. ,  Tj(ti)  =  v(t2)  =  0. 

-1 
Solving  for  C  and  C   by  the  methods  described  in 

Appendices  I  and  II,  v/e  obtain  the  following  expression  for 

u  and  v  evaluated  at  t1    and  t2  respectively: 


•    u(t  )  =  0  =  emi(T~t3>  a-d+r+2kc]  +  em2^"^    [-a-»d°fr-2kc 


(7)  < 


i. 


2r 


-2b-t-k(a-d+r) 


(ta)  =  0  =  ^iU^  |2b±k^a±d±rl  +  em2(T-t2)  : 

A  simplification  of  equations  (7)  can  be  made  by  using  the 
relation  r  =  ma  -  m2 .   The  reduced  form  is  then 

d  +  2kc  -  r 


(8) 


r  pr(T-t4)   _  a. 

a  ~  d  +  2ke  +  r 


r(T-tg)     2b  -  kr  -  k(a-d) 

2b  +  kr  -  k(a-d) 

If  we  now  consider  the  ratio 
sr(T-ti) 


(9) 


_   -r(t2-t!)    a~d+2kc-r  2b+kr-k(a-d) 
r(T~t2)  a~d+2kc+r  '  2b-kr-k(a-d) 


which  can  be  reduced  to 
(10)   er(t2"tl}   =  §-[ 


a  -  d  -  r)  m   b  -  k(a-d)  -  k^c 
a  -  d  +  r)   b  -  k(a-d)  -  k2c  ' 


we  obtain  the  following  relation  between  tA  and  t2 : 


(ID   tg-tj 


log,  (a 


r   '°e  la  -  d  +  r 
that  is,  the  difference  between  the  switching  times  t±   and  t2 
on  a  minimum-time  trajectory  is  an  invariant  cf  the  dynamic 
system. 
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On  the  other  hand,  switching  curves  or  surfaces,  which 
have  received  much  attention,  seem  to  have  no  particular 
significance.   They  are  not  implicit  in  the  differential 
equations,  but  depend  on  the  particular  problem,   lypical 
switching  curves  are  given  in  Illustration  2  for  the  case 
where  the  trajectory  begins  at  the  origin.   In  the  same  figure, 
the  clotted  lines  indicate  the  boundaries  to  the  points  of 
static  stability. 

The  concept  of  using  the  adjoint  system  to  calculate 
differentials  of  end  values  is  apparently  due  to  bliss  in 
his  work  in  ballistics  [£]<>   In  the  problems  solved  here, 
an  explicit  soluticn  to  the  adjoint  system  was  obtained, 
which  made  possible  the  integral  formulas  for  the  end-point  X(T) . 
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AirlEKLIX  I 

Solution  of  a  System  of  Homogeneous  Linear  First-Urder 
Differential  Equations  with  Constant  Coefficients. 

Suppose  we  have  the  differential  equation 

(1)  I   =  BD 

where  B  is  an  n  x  n  matrix  with  constant  elements  b., 
and  v/ith  real,  distinct  eigenvalues,. 

Let  us  assume  a  solution  of  (1)  of  the  form 

(2)  U  =  emt  C 

where  C  is  an  n  x  1  column  matrix  with  constant  elements. 
Differentiating  (2)  with  respect  to  t,  we  get 

(3)  U  =  m  emt  C, 

Substitution  of  (3)  and  (2)  into  (1)  gives  the  following 

result: 

f a\  mt  ^     -c   mt  n 

(4)  me   C  =  B  e   C 

and,  since  e   is  a  scalar  quantity  appearing  en  both 
sides  of  (4) ,  we  have 

(5)  m  C   =  B  Co 

Therefore,  if  m  is  an  eigenvalue  of  matrix  B,  then  C 
must  be  a  corresponding  eigenvector „   Since  there  exist  n 
distinct  eigenvalues  m1?  m2,  0->o»,  mn  of  matrix  B,  there 
are  n  independent  solutions  of  (1)  of  the  form 

(6)  U^   =  em^   C5  (j  =  l,2,..,,n) 

i  T 

where  CJ  =  [ c-,  .  cf. .  0  „  .  c^  .  ]   is  the  eigenvector  corres- 
xj   dj      nj 

ponding  to  eigenvalue  m*,  and  UJ  =  [uj  uj!  <>„<,..  u*?  ]  ° 

lc   In  this  section  only?  ra*  is  associated  v/ith  solutions 
for  Uo 
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The  general  solution  of  (1)  is  then  a  linear  combination 
of  the  bJ,  i.e. , 

(7)    U  =  £k.  U5  =  [  k.e^*  Cj 

where  the  k.  are  arbitrary  constants0 

Componentwise,  (7)  may  be  written  as 

(S)    «±  =  I  kje^*  oi3  (i  =  X,2,...,n). 
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determination  of  the  Components  of  Eigenvector  C  . 

je  have  the  adjoint  system  of  differential  equations 

(1)  u   =   -aT  U 

where  A  is  an  n  x  n  matrix  with  constant  elements  a.,  and 
with  real,  distinct  eigenvalues. 

From  Appendix  I  v/e  knew  that  the  solution  of  (1)  is 
of  the  form 

(2)  U  =  Z  k.  e"121^  Cj 

T         i 
where  the  -m .  are  the  eigenvalues  of  matrix  -A  and  the  Cr 
0 

T 
are  the  corresponding  eigenvectors  of  -A  .   The  k.  are 

j 

arbitrary  constants. 

Differentiating  (2)  with  respect  to  t  we  get 

(3)  u   =   -  7  m.k.  e'™^   C^ 

th 
The  i   component  of  (3)  can  be  written  as 

*  T"        ~m  •  t 

(4)  u.   =   -  I    m.k.  e   3   c. . ,        (i=l,2,...,n). 

th 
vie   can   also  write   the   i        component  of    (1)    as 

\t>)         u  •      =   —   a-j  ^^-j    ~   ap»Up   —    00..0    -   ^riT  ^V '      vi=i>^y»oo5njo 

The  i   component  of  (2)  can  likewise  be  written 
as 

(6)  u-  r  £k.  e"111^  c..,  (i  =  1,2,. ..,n). 

Substituting  (6)  into  (b)  we  get 

/<-7"\  r  V  i   -nut    x       /  T"  i   -m.it    s 

(7)  «!-"«!!(  I  kj«   J  c13)  -  a£i(  £  k.e   Oc^) 

J  J 


•  e  o  o 


Since  (4)  and  (7)  represent  equivalent  expressions  of 
the  u-  we  can  equate  the  corresponding  coefficients  of 
the  k -e  J  .   For  example,  equating  the  corres  ending 

J 

4- 

coefficients  of  kie~  *  we  get 

(&)    -niiCii  =  -aiiCn  -  a2ic2a  -  . ,.  *  anicm  >  (i~  lj2,..n). 

Therefore,  writing  cut  all  n  of  the  equations  obtained  by 

equating  corresponding  coefficients  of  kje"'  1   we  get 

-niiCn  =  -anon  -  a2ic21  -  ....    -   aniCni 

-niic2i  =  -a12Cn  -  a22c21  -  . ...  -  an2cni 
(9) 


~micni  ~  ~amcii  "  a2nc2i  -  <><><»<>  -  anncni 
which  can  be  written  as  the  following  matrix  equation; 


(1C) 


an -mi.    a21 

a12      a2  2~ml 


ani 

"en" 

"  0  " 

an2 

C21 

= 

0 

amTmi 

cm 

! 

I 

0 

am     a2n 

Denote  hy   A  the  determinant  of  coefficients  in  (10). 
For  a  ncntrivial  solution  of  the  cj_i  tc  exist  it  is  necessary 
that  A  =  0,  ioe,,  that 
(11)    I  A1  -  niilj  =  0o 

a  scluticn  of  the  c^i  consists  of  the  ccfactors  of  the 

T 

elements  of  a  row  of  the  matrix  [A  -  mil],  see  Bellman  [  6 ]„ 

This  is  the  same  as  the  cofactors  of  a  column  of  [A  -  cijl], 
That  is,  the  components  of  eigenvector  C1  are  in  the  same 
ratio  as  the  cofactors  of  the  elements  of  a  column  of  the 
characteristic  matrix  [A  -  mal]0 
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Similarly,  the  components  of  the  remaining  eigenvectors 
C  can  be  found  by  using  successively  the  cofactcrs  cf  the 
elements  of  a  column  of  the  matrices  [A  -  m.l]. 
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Illustration  1 
Computer  Flew  Chart  for  Solution  of  N   Order  Minimum -Time  Problem 


Start 


a,   Form  Characteristic  Equation 
of  matrix  A0 


b.   Find  eigenvalues  in.  of  A 

j 

(roots  of  Characteristic  Eq„) 


i     x 

c.   Form  eigenvectors  C   of  -A  o 


CJ  =  [Cii  c 


2    ^j  •••  ^3 


T 


d.   Form  C,  the  matrix  of  eigen- 
vectors of  -Al'o 


e.   Invert  C  to  get  C 


■] 


f.   Load  initial  values  of  the 
input  parameters. 


" 
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a 

!' 
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0 
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\* 

Initial   |f 

Parameters 

^ jbi ,b2  > • • • I 


1 
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SPECIAL  SYMBOLS  AND  NOTATION 
Symbol  Meaning 


T 

X  =  [xi  x2  . .  xl     Matrix  of  coordinates, also  celled 

System  state  function. 
X  =  I  x.e.  Same,  considered  as  a  vector. 


i   i  i 


R  =  xe°!  +  ye2         Same9  considered  as  a  vector. 
F  =  [fi  f2  «o  f  ]     Forcing  (Contrcl)  function. 
F  -  I  f.e.  Same,  considered  as  a  vector. 


I  fiei 


il 


U  =  [ui  u2  oo  u  ]     Solution  to  adjoint  system  of 


n 


differential  equations 


U   =  i   u.e.    =  I  b.    ,U.   Adjoint  vector 
ill        j      j -J.   j 

(b0  =  1) 

UJ=  [ui  ^2  ••  u^]     j "  solution  of  adjoint  system 


n' 


corresponding  to  eigenvalue  -m. 

T    fh 

U.=  [u. j  u2i  ,.u  .]    j   "special"  adjoint  solution 


such  that  u.  .(T)  =  <5  n- ^  =  j'l,  i=j 
ij        13    |  0)  ^ 


Uj  =  I  ui iei  Same,  considered  as  a  vector. 

C  =  [ch  c2n  .  .c  . ]   Eigenvector  of  -A"1"  corresponding 


ou     cj.tcn»diuc     =m  .  o 

T 


C  =  [C1  C2 o . .  C  ]     Matrix  of  eigenvectors  of  -a* 
(U,F)  =  I  u.f.        Scalar  product  of  adjoint  vector 


ill 


and  forcing  vector 0 


e1,e2,...,e  Base  vectors  of  Cartesian  coordinate 

set  in  n= space. 

A  =  (a-, )  n  x  n  matrix  of  coefficients. 

B  =  (b-O  n  x  n  matrix  of  coefficients. 

In. 

[a  -  ml]  Characteristic  matrix  of  a. 

[ii  -  ra-I]  Characteristic  matrix  of  A 

associated  with  eigenvalue  m.. 

| A  -  ml |  =0         Characteristic  equation  cf  A, 

5  3 


Symbol 


Meaning 


th  . 


Diagonal  matrix  with  i  diagonal 

-,     ,   mi(T-t) 
element  e  J 


X(0)  =  (xi,x2r,^)   Initial  state  of  dynamic  system 


n'o 


X^  =  (X-, -,Xg«, . .  ,X  «)  Terminal  state  of  oynamic  system 

(specified) 

X(t)  =  (Xi ,x2, . o ,x  )n.  Terminal  state  of  dynamic  system 

at  time  t  =  To 


T 
f 


Terminal  value  of  time;  also 
transpose  when  used  with  matrix. 

Subscript  f  denotes  final  value. 


! 

c .  • 
AX, 


kif 


x.(T) 


L_  =  (U.F  ),  t,<t<t 


Bar  over  a  symbol  indicates  that 
it  is  a  vectoro 


Element  of  C   . 

Error  in  terminal  value  of  i 
component . 


th 


k+1 


*Xi  =  «£«  ♦  gj*bl  +  ....  ♦&*  «Vi 


t   ]t 


ik 


jump  discontinuity  in  variable 
at  time  t., 

IK 


£fl]t«,.  =  W.^  -   ^(tik-O) 


(5j.[F]) 

'ik 


'ik 


'  T   t^  +0 
J    ik 


(V\k-o 


iVfi] 


is.  =  L^ie      k2e  ^ 


k   switching  time  of  f. 
k  e°mntlT 


'ik 


£4 
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